INFORMATION  THEORETIC  RADAR  WAVEFORM  DESIGN  FOR  MULTIPLE  TARGETS 


Amir  Leshem,  Oshri  Naparstek 

School  of  Engineering, 
Bar-Ilan  University, 
Ramat-Gan,  Israel 


ABSTRACT 

In  this  paper  we  describe  the  optimization  of  an  information 
theoretic  criterion  for  radar  waveform  design.  The  method  is  used  to 
design  radar  waveforms  suitable  for  simultaneously  estimating  and 
tracking  parameters  of  multiple  targets.  Our  approach  generalizes 
the  information  theoretic  water-filling  approach  of  Bell.  The  paper 
has  two  main  contributions.  First,  a  new  information  theoretic 
design  criterion  for  designing  multiple  waveforms  under  a  joint 
power  constraint  when  beamforming  is  used  both  at  transmitter  and 
receiver.  Then  we  provide  a  highly  efficient  algorithm  for  optimiz¬ 
ing  the  transmitted  waveforms,  by  approximating  the  information 
theoretic  cost  function.  We  show  that  using  Lagrange  relaxation  the 
optimization  problem  can  be  decoupled  into  a  parallel  set  of  low- 
dimensional  search  problems  at  each  frequency,  with  dimension 
defined  by  the  number  of  targets  instead  of  the  number  of  frequency 
bands  used. 
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I.  INTRODUCTION 

The  problem  of  radar  waveform  design  is  of  fundamental  impor¬ 
tance  in  designing  state-of-the-art  radar  systems.  The  possibility  to 
vary  the  transmitted  signal  on  a  pulse-by-pulse  basis  opens  the  door 
to  great  enhancement  in  estimation  and  detection  capability  as  well 
as  improved  robustness  to  jamming.  Furthermore  modern  radars 
can  detect  and  track  multiple  targets  simultaneously.  Therefore, 
designing  the  transmitted  waveforms  for  detecting  and  estimating 
multiple  targets  becomes  a  critical  issue  in  radar  waveform  design. 

Most  of  existing  waveform  design  literature  deals  with  designs 
for  a  single  target.  One  of  the  important  tools  in  such  designs  is 
the  use  of  information  theoretic  techniques,  see  [1]  for  a  review 
of  early  results  of  Woodward  and  others.  Bell  [1]  was  also  the 
first  to  propose  using  the  mutual  information  between  a  random 
extended  target  and  the  received  signal.  His  optimization  led  to 
a  water-filling  type  strategy.  In  his  paper  he  assumed  that  the 
radar  signature  is  a  realization  of  random  Gaussian  process  with  a 
known  power  spectral  density  (PSD).  However,  when  considering 
real-time  signal  design  we  can  use  his  approach  to  enhance  the 
next  transmitted  waveform  based  on  the  a  priori  known  signature. 
Whereas  waveform  design  literature  concentrated  on  the  estimation 
of  a  single  target,  modern  radars  treat  multiple  targets.  Therefore, 
the  development  of  design  techniques  for  multiple  targets  is  of 
critical  importance  to  modern  radar  waveform  design. 

Recently  a  great  interest  has  emerged  in  MIMO  radars,  where 
multiple  transmit  and  receive  antennas  are  used  with  large  spatial 
aperture  to  overcome  target  fading  ,  see  [2]  and  the  references 
therein  for  a  review  of  statistical  MIMO  radar  literature.  Much  less 

The  work  of  A.  Nehorai  was  supported  by  the  Department  of  Defense 
through  the  Air  Force  Office  of  Scientific  Research  MURI  grant  FA9550- 
05-1-0443  and  AFOSAR  Grant  FS9550-05-1-0018.  Contact:  Amir  Leshem, 
leshema  @  eng .  biu  .ac.il 


Arye  Nehorai 

Department  of  ESE, 
Washington  University, 
St.  Louis,  Missouri 
U.S.A 


has  been  done  on  MIMO  waveform  design.  Two  papers  related 
specifically  to  waveform  design  in  the  MIMO  context  are  [2], 
[3].  Yang  and  Blum  applied  MIMO  point-to-point  communication 
theory  to  design  radar  waveforms  by  water-filling  the  power  over 
the  spatial  modes  of  the  overall  radar  scene  (channel).  They  also 
showed  that  optimizing  the  non-causal  MMSE  and  optimizing 
the  mutual  information  leads  to  identical  results.  This  finding 
provides  another  justification  for  using  the  maximum  mutual  in¬ 
formation  criterion  for  the  radar  waveform  design  problem.  Their 
work  is  a  novel  extension  of  the  work  of  [1].  However,  one 
should  note  that  by  water-filling  with  respect  to  the  spatial  modes, 
higher  power  is  allocated  to  the  stronger  targets.  This  approach 
is  reasonable  when  using  a  single  target  through  several  remotely 
located  antennas.  However,  this  approach  is  not  always  desirable, 
when  tracking  multiple  targets.  De  Maio  and  Lops  [3]  proposed 
design  criterion  for  space  time  codes  for  MIMO  radars  based  on 
mutual  information.  They  also  analyzed  the  detection  probability 
of  these  techniques  under  the  statistical  MIMO  diversity  model  that 
assumes  independent  scattering  towards  each  of  the  MIMO  systems 
component  and  point  targets. 

The  approach  proposed  in  this  paper  is  different.  We  are  inter¬ 
ested  in  reception  and  transmission  for  tracking  multiple  extended 
targets,  by  using  the  insights  provided  by  multi-user  information 
theory  instead  of  the  point-to-point  MIMO  approach.  These  insights 
are  applied  here  for  the  context  of  coherent  phased  array  receivers 
that  are  capable  of  transmitting  independent  signals  simultaneously, 
as  well  as  for  optimizing  the  waveforms  for  extended  targets.  We 
assume  high  range  resolution  and  that  the  various  extended  targets 
are  treated  as  independent  signals  that  need  to  be  estimated.  In  the 
optimization  process  we  provide  priorities  through  a  set  of  priority 
vectors.  A  linear  combination  of  the  mutual  information  between 
each  radar  beam  and  its  respective  target  is  optimized.  This  leads  to 
a  highly  complicated  optimization  problem.  However,  by  assuming 
linear  pre-  and  post-processing  and  an  independent  estimation  of 
the  targets,  we  are  able  to  reduce  the  waveform  design  problem 
to  a  problem  similar  to  that  of  the  centralized  dynamic  spectrum 
allocation  in  communication.  Furthermore,  recent  advances  in  con¬ 
vex  optimization  (see  [4]  and  the  references  therein)  open  the  way 
to  design  techniques  specifically  tailored  for  radar  waveforms  that 
would  be  suitable  for  estimating  the  parameters  of  multiple  targets. 
While  we  concentrate  on  the  phased  array  coherent  reception, 
statistical  MIMO  modeling  might  be  incorporated  into  this  context 
where  each  transmitter  optimizes  its  transmit  spectrum  towards  the 
various  independent  targets. 

The  paper  has  two  main  contributions:  First  we  extend  Bell’s  re¬ 
sults  to  the  design  of  multiple  transmit  waveforms,  each  optimized 
towards  a  specific  target  where  the  transmitter  employs  multiple 
beamformers  as  well  as  receiver.  Finally,  an  optimization  algorithm 
is  proposed.  We  show  that  using  duality  theory  the  problem  can 
be  reduced  to  a  search  over  a  single  parameter  and  parallel  low¬ 
dimensional  optimization  problems  at  each  frequency.  Interestingly 
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even  though  the  proposed  design  criterion  for  multiple  waveforms 
is  non-convex,  strong  duality  [4]  still  holds,  which  allows  us  to 
solve  the  simpler  dual  problem. 

II.  TARGETS  MODEL 

In  this  section  we  describe  the  extended  targets  model.  While 
classical  radar  target  models  assume  far-held  point  source  targets. 
This  is  indeed  the  case  when  the  radar  pulse  is  relatively  narrow- 
band  so  that  the  range  span  of  the  target  is  well  within  a  single 
range  cell.  In  contrast  to  these  point  source  models,  many  modern 
radars  are  often  capable  of  transmitting  very  wide-band  pulses 
or  alternatively  use  very  wide-band  compressed  signals.  In  this 
case  delays  across  the  target  are  similar  in  nature  to  multipath 
propagation.  This  results  in  a  complex  target  impulse  response. 
Some  examples  of  wide  band  responses  of  airplanes  and  missiles 
can  be  found,  e.g.,  in  [5].  Under  these  conditions  the  targets  are 
called  extended  targets,  which  are  the  focus  of  the  current  paper. 
Models  for  such  targets  have  been  used  e.g.,  in  [1].  Extended  targets 
naturally  appear  in  imaging  and  high  range  resolution  applications 
[6]  where  the  radar  signal  bandwidth  is  sufficiently  large  so  that  the 
target  is  not  contained  in  a  single  range  cell.  Such  target  models 
were  already  described  by  Van  Trees  [7]  where  they  are  termed 
range  selective  targets.  Extended  targets  typically  have  multiple 
reflection  centers,  each  with  independent  statistical  behavior.  The 
target  impulse  response  (TIR)  is  therefore  modeled  as 

L 

n(t)  =  -  re)  (1) 

£=1 

where  ti  <  2 d/c,  d  is  the  radial  span  of  the  target  and  c  is  the 
speed  of  light,  m  are  the  individual  random  reflection  coefficients. 
These  coefficients  can  be  modeled  either  deterministically  or  using 
the  extended  Swerling  x2  models  [8].  The  temporal  variability  of 
the  target  response  is  mainly  determined  by  the  speed  of  the  target 
and  the  carrier  frequency.  The  reflected  radar  signal  is  given  by 

rT  max 

y(t)  =  /  K(r)s(t  -  t )dr  (2) 

Jo 

where  k(t)  is  the  TIR  and  s(t)  is  the  radar  signal.  Since  the  targets 
have  non-trivial  impulse  response,  we  can  consider  also  the  target 
frequency  response  (TFR)  given  by 

h(f)  =  /  K(T)e-*"frdT 
Jo 

For  stochastic  target  models,  we  will  be  interested  in  the  PSD  of  the 
TIR  which  now  becomes  a  stochastic  process.  Significant  amplitude 
deviations  will  only  appear  for  extended  targets.  Typically  we  will 
sample  the  frequency  domain  and  assume  that  h(f)  is  given  at  a 
set  of  K  equally  spaced  frequencies  h(fk)  :  k  =  1, ...,  K. 

III.  PHASED  ARRAYS  TRANSMIT  AND  RECEIVE 
BEAMFORMING 

Our  main  interest  in  this  paper  is  with  phased  arrays  which 
use  both  transmit  and  receive  beamforming.  We  now  provide 
the  basic  model  for  transmit  and  receive  beamforming  for  mul¬ 
tiple  targets.  We  can  typically  assume  that  the  array  manifold 
is  independent  of  frequency.  This  holds  as  long  as  the  transmit 
signal  bandwidth  is  small  relative  to  the  carrier  frequency.  In 
this  section  we  will  maintain  this  assumption,  but  we  will  use 
the  more  general  formulation  in  the  following  sections.  Assume 
that  we  have  a  phased  array  radar  capable  of  transmitting  and 
receiving  simultaneously  L  beams.  Each  beam  is  characterized 
by  transmit  beamforming  vectors  (um,m  =  1,...,L)  and  receive 
vectors  (w m,ra  =  1,...,L).  The  baseband  signals  Si(t)  that  are 
transmitted  over  the  respective  beams  are  multiplied  by  the  transmit 


beamforming  vectors  and  linearly  combined  to  form  the  baseband 
transmit  vector 

L 

t  (t)  ^  UmSm(t)  (3) 

m=  1 

Let  a(0)  be  the  array  manifold  of  the  array  towards  direction  9. 
The  transmitted  signal  is  reflected  at  a  target  with  direction  0  and 
range  R  is  given  by 

y(0,  *)  =  -]=;/  K(T)a(0)*t  dr  (4) 

where  we  neglect  the  free  space  attenuation  across  the  target  (since 
CTmax  <<  R)-  The  reflection  of  a  target  at  direction  6  is  received 
by  the  array  as  x(t)  =  1/Ra.(0)y(0,t).  Assuming  that  we  have  L 
targets  with  directions  (0i,...,0l)  and  ranges  (Ri  :  £  =  1, ...,  L) 
we  obtain  that  the  received  signal  is  given  by 

x(i)  = 

Etx  EU 1  ^  Kl(T)sm  (t-(%+  t))  dr 

(5) 

where  =  a(^)a(^)*.  To  enhance  the  signal  to  noise  ratio  by 
suppressing  directional  interference  and  other  targets  side-lobes  we 
apply  L  transmit  beamforming  vectors  (w^  :  £  =  1,  ...,L)  to  the 
received  signal  resulting  in 

ze{t)  =  w?x(t)  (6) 

This  is  the  standard  way  to  decouple  the  estimation  between 
azimuth  cell,  since  it  greatly  reduces  the  number  of  targets  that  need 
to  be  estimated  jointly.  Using  (5)  and  translating  to  the  frequency 
domain  we  now  obtain 

L  L  1 

ze(f)  =  '22'22  #2  M/)  (w«a(6><0)  (a(0c)*um)  sm(f)  (7) 

£=1  m=l  ^ 

where  he(f)  is  the  Uth  target  frequency  response.  To  simplify 
notation  from  this  point  on  we  will  assume  that  is  included 

in  the  target  signature.  When  the  targets  are  resolved  in  range 
or  in  angle  we  can  separate  them  in  the  time  domain  or  using 
receive  beamforming,  which  means  that  only  certain  range  cells 
will  include  target  information.  This  will  imply  that  each  ze(t)  is 
subject  to  only  receiver  and  clutter  noise.  When  targets  are  partially 
overlapping  both  in  range  and  angle  (see  e.g.,  Gini  et.al  [9])  each 
beam  contains  residual  interference  from  other  targets.  In  this  case 
the  noise  PSD  contains  contributions  from  other  targets.  The  next 
step  is  a  correlation  of  each  ze(t)  with  se(t)  to  obtain  the  target 
impulse  responses.  These  impulse  responses  can  be  used  to  enhance 
the  transmitted  signal  in  the  next  pulse.  This  can  be  done  by  using 
the  targets  PSD  when  the  target  reflection  centers  (and  therefore 
the  target  signature  PSD)  exhibit  pulse  to  pulse  variations  as  in 
the  Swerling  type  II  models  or  by  using  the  latest  estimate  when 
the  variations  are  sufficiently  small.  The  exact  choice  of  the  model 
depends  on  the  target  velocity,  radar  carrier  frequency  and  PRI  or 
compressed  pulse  duration. 

Since  the  targets  are  selective  in  range,  we  also  obtain  that  certain 
frequencies  are  more  reflective.  This  implies  that  concentrating  the 
transmitted  power  according  to  the  target  frequency  response  is 
beneficial  in  terms  of  the  information  we  obtain  regarding  the  target 
signature. 

III-A.  Multi-target  tracking 

Finally  we  discuss  the  tracking  model,  and  its  relationship  to  the 
signal  design  problem.  In  general  multi-target  tracking  is  a  well 
established  topic  [10].  Our  paper  is  not  focused  on  the  tracking 
itself  but  rather  on  the  adaptive  design  of  the  transmitted  waveform, 
based  on  the  target  parameters.  Therefore  the  design  will  be  affected 
by  the  following  parameters: 


1.  The  azimuth  and  range  cells  that  include  each  target.  These 
influence  the  transmit  and  receive  beamforming  vectors. 

2.  Target  motion  during  the  time  interval  between  pulses  relative 
to  the  carrier  wavelength.  This  parameter  decides  the  statis¬ 
tical  model  of  choice  for  estimating  the  TIR.  If  the  motion 
is  large  compared  to  the  wavelength  then  we  can  use  only 
target  PSD  as  in  the  Swerling  type  II  or  IV,  while  if  the 
motion  is  small  so  that  the  local  reflection  environment  can 
be  considered  static  we  can  use  the  previous  estimate  of  the 
TIR  as  a  predictor  for  the  next  realization. 

Since  our  main  interest  is  in  adaptive  design  of  the  pulse,  we 
shall  assume  a  given  estimate  for  these  parameters,  assume  that 
the  transmit  and  receive  beamforming  vectors  for  each  beam  are 
provided  by  the  tracking  system,  and  limit  our  interest  to  the  radar 
signal  design  problem.  This  is  a  reasonable  approach  since  the 
described  parameters  are  provided  by  existing  systems.  We  will 
also  assume  that  the  radar  control  provides  us  priorities  with  respect 
to  the  various  targets  to  be  tracked.  These  priorities  are  given  by 
a  vector  of  constants.  The  choice  of  these  constant  is  important. 
However,  the  relative  priorities  can  be  determined  from  the  overall 
SNR  estimate  of  each  target  as  well  as  its  temporal  variability, 
which  depends  on  the  target  speed.  Typically  we  would  like  to 
allocate  higher  priority  to  rapidly  moving  targets  or  weak  targets 
that  are  harder  to  track. 

IV.  INFORMATION  THEORETIC  APPROACH  TO 
WAVEFORM  DESIGN 

In  this  section  we  extend  the  waveform  design  paradigm  of  Bell 
[1]  to  the  case  of  multiple  radar  transmitters  and  receivers.  In  order 
to  study  the  trade-off  between  various  radar  receivers,  we  use  a 
linear  convex  combination  of  the  mutual  information  between  the 
targets  and  the  received  signal  at  each  receiver  beam  oriented  at 
that  specific  target. 

We  begin  by  revising  the  received  signal  model.  Assume  that  an 
array  with  p  elements  simultaneously  transmits  L  waveforms.  The 
transmitted  signal  at  frequency  fk  is  given  by 


where  v'i(fk)  =  w}(fk)v(fk)  is  the  received  noise  and  clutter 
component  of  the  Cth  beam.  Let  cr^,(fk)  =  E  \i/'£(fk)\2  Af 
be  the  Cth  beam  noise  power  at  frequency  fk .  After  algebraic 
manipulations  we  can  show  that  the  mutual  information  between 
the  Cth  beam  and  the  tlh  target  at  frequency  fk  is  now  given  by 

'* = 108  (‘ +  A/' 

(13) 

where  the  signal  reflected  from  the  Eth  target  is  denoted  by 

L 

ze(fk)  =  D  w^(/fc)R<(/fc)um(/fc)sm(/fe).  (14) 

m=  1 

while  the  noise  and  inter- target  interference  component  at  the  Eth 
beam  is  given  by 

L 

Ze(fk)  =  EE  WK/fc)Rn(/fc)U»n(/fe)Sm(/fc)+^(/fc)-  (15) 

n^l  m=l 

We  assume  that  the  radar  allocates  one  beam  towards  each  target, 
since  non-linear  joint  processing  of  all  the  beams  would  lead  to  an 
infeasible  receiver.  Therefore,  the  total  mutual  information  between 
the  ^’th  beam  and  the  ^’th  target  is  given  by: 

K 

I  (he;  ze\s)  =  ]T  /,.  (he(fk);  ze(fk)\s(fk))  A /  (16) 

k=  1 

where  z£  =  (zi(fk)  :  k  =  1, ...,  K)  and  h£  = 
{hi(fk)  ■  k  =  1, ...,  K)  are  the  received  signals  using  the 
t?’th  received  beam  and  the  tth  target  signature,  respectively, 
s m  =  [sm(/i),-)sm(/x)]T  are  the  signal  waveform  samples 
directed  towards  the  ra’th  target, 

S  =  [si,...,sl] 


L 

t(fk)  =  J2Mfk)se(fk),  k=l,...,K  (8) 

1  =  1 

where  u e(fk)  are  the  beamformer  coefficients  for  the  E th  wave¬ 
form  designed  for  the  Eth  target  at  frequency  fk,  and  se(k)  is 
the  corresponding  waveform  at  frequency  fk.  We  assume  channel 
reciprocity;  i.e.,  if  the  receive  steering  vector  is  a (0£,  fk),  then  the 
transmitted  signal  arrives  at  the  target  with  channels  a *  (#£,/&)• 
The  signal  reflected  from  the  Eth  target  having  signature  h£  = 
(ht(fk),k  =  1, ...,  K)  is  therefore  given  by 

L 

y e(fk)  =  L  (a* (Oe,  fk)um(fk))  he(fk)smUk)  (9) 

m= 1 

for  k  =  1 AT  (note  that  we  have  used  index  m  to  enumerate 
the  transmitted  waveforms,  m  =  1,...,L,  since  i  is  reserved  for 
the  target).  Hence,  the  received  signal  at  the  array  is  given  by 

L 

x(/fc)  =  ^2  R(fk)um(fk)sm(fk)  +  v(fk),  (10) 

m=  1 


is  the  complete  spatio-temporal  waveform  matrix,  and  s  =  vec(S). 
Assuming  that  the  beamforming  vectors  are  known  the  multiple 
waveform  design  problem  is  now  given  by 

maxs  Yj^=i  (h£;z^|s) 

subject  to  ttl  \S^k\2  <  Tmaxj  1 

where  a  =  [ai , ...,  oll]T  is  the  target  priority  vector.  This  problem 
is  highly  non-linear  in  the  complex  waveforms  S.  Furthermore,  it 
involves  cross-correlations  between  the  waveforms,  and  therefore 
phase  information  plays  an  important  role.  Hence  we  need  to 
design  not  only  the  waveform  spectrum,  but  the  complete  complex 
envelope.  The  dependence  on  the  phase  will  have  a  secondary 
drawback,  since  we  will  not  be  able  to  reduce  the  peak  to  average 
of  the  overall  transmitted  waveform  by  properly  choosing  the 
waveform  phase.  However,  we  will  show  that  in  the  typical  scenario 
of  multiple  beams  in  a  large  phased  array  this  problem  can  be 
approximated  by  a  simpler  spectrum  design  problem. 

V.  WAVEFORM  OPTIMIZATION  FOR  MULTIPLE 
TARGETS 


where  R (/*)  =  Xf=i  R«(/fc)-  and 

R e(fk)  =  he(fk)a(9e,  fk)a* (Be,  fk).  (11) 

Assume  that  a  beamformer  w e(fk)  is  used  to  receive  the  Eth 
target,  resulting  in 

Zi(fk)  =  w|(/fc)]T^=iR(/fc)um(/fe)Sm(/fc)  +  v'e(fk) 


Using  certain  approximation  of  the  mutual  information  and 
relying  on  the  properties  of  the  transmit  and  receive  beamformers 
we  can  show  that  the  mutual  information  (16)  can  be  approximated 
by 


K 

i{  h«;zc|p)  =  y>g 

k= 1 


_ P£,k  _ 

E m^£  9£,m(fk)Pm,k  +  CT^/  (fk) 


(18) 


where  pm  =  [pm,i,  ••• ,Pm,K]T  is  the  power  allocation  for  the  ra’th 
target, 

p  =  [Pl,...,pL] 

is  the  total  power  allocation  matrix,  and  p  =  vec(P).  The  constants 
gi^rn  are  defined  by 

9e,m(fk)  =  wJ(/fe)Rm(/fe)um(/fe) 

Pm,k  —  \sm(fk)\  A/, 

and  include  all  the  prior  information  regarding  the  target  signatures 
and  the  channels. 

The  problem  (17)  can  now  be  simplified  to 

maxP  XAi  (h<;  z^|p)  (19) 

subject  to  Y^i=1  Yjk=iVi,k  <  Tnax 

To  solve  the  multiple  waveform  design  problem,  we  should  note 
that  (19)  is  a  generalized  (non-convex)  monotropic  optimization 
problem,  since  the  summands  of  (18)  are  not  concave  functions. 
While  this  is  a  non-convex  optimization  problem  we  will  show 
how  it  can  be  solved  efficiently  using  duality  theory. 

Applying  duality  theory  we  obtain  that  the  Lagrangian  dual 
function  is  now  given  by 

K 

Ld( A)  =  infp  —  y  ^  Lk  (p*.,  A)  A /  +  A Pmax  (20) 

k=i 

where 

L 

Lk(pk ,  A)  =  ^  OLfJk  (M/fc);z*(/fc)|pfc)  A/  +  AlTpfc.  (21) 
£=1 

The  dual  problem  now  becomes 

maX  V  ifc(Pfc>  “  ^-Pmax  I  •  (22) 

A-°  \k=lPk  ) 

Note  that  unlike  the  case  of  a  single  waveform,  we  will  have 
multi-dimensional  parallel  optimization  problems.  However,  this 
problem  has  two  significant  simplifications:  The  dimension  L  of 
each  problem  is  much  smaller  than  the  typical  number  of  frequency 
bins.  Second,  the  problem  is  unconstrained,  which  is  a  major 
simplification  in  the  non-convex  problem. 

We  can  now  solve  (22)  using  bi- section  search  for  A  and  solving 
the  parallel  problems  at  each  frequency  given  any  specific  value  of 
A.  This  is  done  using  standard  unconstrained  optimization  tools. 
While  the  complexity  is  still  large,  it  is  still  linear  in  the  number 
of  frequency  bins.  Furthermore  our  functions  are  smooth,  and  the 
gradient  and  Hessian  are  rational  functions.  This  can  be  exploited 
in  solving 

pk  =  inf  Lfc(pfc,  A) 

P  k 

In  simulations  we  will  show  how  these  problems  can  be  solved 
efficiently  when  the  total  number  of  variables  is  much  above  100. 

VI.  SIMULATIONS 

In  this  section  we  show  some  simulated  results,  we  assume  that 
two  waveforms  are  transmitted  by  an  omni-directional  equispaced 
linear  phased  array  with  10  elements  (A  spacing)  and  received  by 
the  same  array.  The  target  directions  in  our  simulation  were  70° 
and  80°.  Waveforms  bandwidth  was  80MHz.  We  have  used  100 
frequency  bins,  corresponding  to  800  kHz  frequency  resolution. 
The  receive  beamformers  used  were  based  on  MVDR  beamformer, 
and  the  transmit  beamformers  were  classical  beamformers  directed 
towards  the  targets.  Target  signatures  were  Gaussians  corresponding 
to  targets  of  length  17m  and  10m  respectively  as  shown  in  Figure 
1.  The  priorities  used  for  the  experiment  were  a\  —  0.4,a2  = 


Targets  impulse  response 


Target  Frequency  response  (normalized) 


Fig.  1.  Gaussian  modeled  target  responces  for  17m  and  10m  long 
targets. (a)  Target  impulse  response,  (b)  Target  frequency  response. 


0.6  with  the  higher  priority  given  to  the  weaker  target,  the  two 
targets  were  chosen  to  be  spatially  separated.  Figure  2  shows  the 
transmitted  PSD  for  both  targets.  In  this  example  we  can  see  that 
the  algorithm  transmits  for  both  targets  with  a  very  large  frequency 
overlap.  This  could  be  explained  by  the  received  PSD  shown  in 
Figure  3. The  targets  do  not  strongly  interference  with  each  other 
and  the  interference  spectrum  is  relatively  flat  over  all  frequencies. 
This  implies  that  the  algorithm  can  transmit  for  both  targets  in  the 
frequencies  where  their  SINR  is  high  even  if  the  frequencies  are 
overlapping,  without  losing  information  due  to  interference  between 
the  targets. 


Desinged  waveforms  for  two  targets 


Fig.  2.  Transmitted  waveforms  towards  the  targets. 

Finally  we  have  studied  the  information  region  of  the  two  targets 
and  compared  to  the  case  where  no  spectral  shaping  is  applied  to 
the  transmitted  pulse.  The  results  are  presented  in  figure  4.  Using 
flat  spectrum  causes  a  loss  of  100%  for  the  weak  target  compared 
to  the  case  where  the  design  is  according  to  the  weak  target  profile. 
However  choosing  a  —  0.5  leads  to  performance  enhancement  of 
33%  for  both  targets. 

VII.  CONCLUSIONS 

In  this  paper  we  discuss  the  optimization  of  multiple  waveforms 
for  multiple  targets  under  joint  power  constraint.  This  type  of 
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Fig.  3.  Received  waveforms,  after  beamforming. 


Fig.  4.  Information  region,  for  two  extended  targets 


waveform  design  is  suitable  for  unresolved  extended  targets.  We 
have  derived  computationally  efficient  algorithm  and  presented  the 
result  of  the  optimization  in  simulations.  Further  results  as  well  as 
design  of  a  single  waveform  optimized  for  multiple  targets  can  be 
found  in  [11]. 
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